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Abstract. We analyse the ultraviolet variability time scales in a sample of 15 Type 1 Active Galactic Nuclei (AGN) observed 
by IUE. Using a structure function analysis, we demonstrate the existence in most objects of a maximum variability time 
scale of the order of 0.02-1.00 year. We do not find any significant dependence of these maximum variability time scales 
on the wavelength, but we observe a weak correlation with the average luminosity of the objects. We also observe in several 
objects the existence of long-term variability, which seems decoupled from the short-term one. We interpret the existence of a 
maximum variability time scale as a possible evidence that the light curves of Type 1 AGN are the result of the superimposition 
of independent events. In the framework of the so-called discrete-event model, we study the event energy and event rate as a 
function of the object properties. We confront our results to predictions from existing models based on discrete events. We show 
that models based on a fixed event energy, like supernova explosions, can be ruled out. In their present form, models based on 
magnetic blobs are also unable to account for the observed relations. Stellar collision models, while not completely satisfactory, 
cannot be excluded. 
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1. Introduction 

The UV excess of the spectral energy distribution of Type 
1 Active Galactic Nuclei (AGN), the blue bump, reflects the 
fact that a very large fraction of the energy is released in th e 
wavelength domain ~ 300 to 5600 A (see e.g. lKrolikl (H999). 
Fig. 7.10). Conventional accretion disk models are able to ac- 
coun t satisfactorily for the r ough shape of the blue bump (but 
see iKoratkar & Blaesl J 19991)1. but th ey fail to explain the vari- 
ability properties ( Courvoisier & Clavell99lh . a key to the un- 
derstanding of the AGN phenomenon. 

Th is difficulty led several a uthors ICid Fernandes et ail 
1 19961). IPaltanj & CourvoisieJ 1 19971) (hereafter PC97), 
lAretxaea et alJ 1 19971) and |Cid Fernandes et alJ (2000) (here- 
after CSV00); see also Uretxaga & Terlevichl ( 1 1994 ) to 
consider a more phenomenological approach based on the 
discrete-event model, which provides a simple explanation 
for the variability: The variability is the result of a su- 
perimposition of independent events occurring at random 
epochs at a given rate. The motivations for the discrete-event 
model are twofold: Temporal analysis allows to constrain the 
event properties, while its generality leaves room for a large 
variety of physical events. It can be a reasonable approxima- 
tion for mode l s like starburst lAretxaga & Terlev ichl 1 19941 
lAretxaga et alJ Il997l) . stellar collisions ICourvoisier et alJ 
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Il996t iTorricelli-Ciamponi et alJ l2000h. o r magnetic blobs 
above an accretion disk lHaardt et alll994l) . 

In this paper, we use data from the International Ultraviolet 
Explorer (IUE) covering about 17 years to determine the UV 
characteristic time scales in a sample of Seyfert 1 galax- 
ies and QSOs, using a m ethodology similar to that used in 
ICollier & Peterson! 1200 lb (hereafter CP01). While CP01 con- 
centrated on the measure of time scales shorter than 100 days 
by selecting short portions of the light curves in which the time 
sampling was denser, we use here the full available light curves 
(on average 16.5 years), highlighting a wider range of time 
scales, and extend the sample to 15 objects. Furthermore, using 
12 wavelength windows between 1300 and 3000 A, we investi- 
gate for the first time the existence of a wavelength dependence 
of the variability time scale. 

We interpret the variability properties of our objects in 
terms of discrete-event model, and we study their parameters 
as a function of the object properties. Our approach is simi- 
lar to that of CSV00, although we use data from IUE gath- 
ered during almost 17 years, while they use optical data cov- 
ering about seven years. Furthermore, the sample of CSV00 
is composed of PG quasars (median z: 0.16) while ours is 
mainly composed of Seyfert 1 galaxies at a much smaller red- 
shift (median z: 0.033). Their observations thus not only cover 
about three times less time than ours, but the observation dura- 
tions are further diminished in the observer's fra me. Finally, as 
variability increases towards shorter wavelengths ( Kinney et al. 
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199 it IPaltani & Courvoisierl Il994l Idi Clemente etail Il996l 
Trevese & Vag netti 2002), the study of the variability is more 
efficient in the UV than in the optical. 

2. The concept of discrete-event model 

The forma lism of the discrete-event model was mainl y de- 
velope d in ICid Fernandes et all lll996h . PC97, lAretxaga et alJ 
ill 9971) . and CS V00. In the discrete-event model, the variability 
is due to the superimposition of independent events, occurring 
at random epochs, on top of a possible constant source. In the 
simplest form that we use here, all events are identical. 

The total luminosity density at wavelength A can be ex- 
pressed as: 



L A (t) = Yu ei(f ~^ + Ci 



(1) 



where ex(t — f,) is the light curve of event i, initiated at f,, and 
Cx reflects the possible contribution of a steady component. 
The distribution of f,- is assumed to be Poissonian. Using N for 
the event rate and E A = j e A (t - f,)df for the energy density 
released by the 2 th event, the average luminosity reads: 



NE X + Cx 



(2) 



Parameterizing the event using its duration 2/ix and its ampli- 
tude at maximum Hx, we have: 



L A = Nk L H A 2fi A + C A 



(3) 



where &l is a constant depending on the event shape. The 
variance of Lx(t) was calculated by PC97 (see their equation 
(A5)), which reduces to Var(L,i) oc N. Including Hx and 2/ix in 
Eq. (B2) of their paper, we find: 



Var(L,0 = Nk y H 2 A 2nx 



(4) 



where ky is a constant depending on the event shape. From Eqs. 
and ©, the event amplitude at wavelength A reads: 



H A = 



k L Var(L i ) 
kvTx-Cx 



(5) 



In Sect. 0] we describe a method to estimate 2fix- As al- 
ready noted by CSV00, our system will not be closed, as only 
three parameters can be measured: (Lx, Vax(Lx), 2/j.x), while 
our model requires the knowledge of six parameters (N, Ex, 
2fix, Cx, kh, and ky). Fixing the event shape determines k^ and 
ky. Cx is unknown, but constrained in the range < Cx < 
min, Lx(t), see Sect. 15. 21 The system can therefore be solved for 
these limiting cases. Under these assumptions, we can therefore 
derive the energies and rates of the events from the light curves. 

3. The IUE light curves 

3.1. Data selection 

We selected all the Type 1 AGN spectra available in 
early December 2001 in the INES (IUE Newly Extracted 
Spectra) v3.0 database at VILSPA/LAEFF, which used 
a new noise model and background determination 
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Fig. 1. Average IUE spectrum of the Seyfert 1 galaxy Mrk 335. 
The position of the 12 spectral windows is indicated by the gray 



jRodrfguez-Pascual et a Il999h . We extracted the objects 
monitored for several years for which at least 20 large aper- 
ture, small dispersion observations have been performed with 
the SWP instrument (1150-1950 A). These conditions were 
imposed by the temporal analysis (see Sect. 14. 1> . Table^gives 
a list of the selected objects with their common names and 
redshifts. We finally have in our sample 13 Seyfert 1 galaxies, 
one broad-line radio galaxy (BLRG), and one quasar. 

When building the light curves, FITS headers of all spec- 
tra were carefully checked for anomalies. We excluded spectra 
which were affected by an objective technical problem stated 
in the FITS headers (e.g., no significant flux detected, object 
out of aperture, no guiding, no tracking). Spectra for which the 
pointing direction was farther than 10" from the object posi- 
tion were also discarded. For 3 C 273, the selecte d spectra cor- 
respond to the list described in Tiirler et al. ( 1999). 

3.2. Light curves 

For each object, we built 12 light curves in 50 A spectral 
windows starting at 1300 A, 1450 A, 1700 A, 1950 A, 2100 A, 
2200 A, 2300 A, 2425 A, 2550 A, 2700 A, 2875 A, and 2975 A 
in the rest frame of the object, avoiding contamination by 
strong emission lines. Fig. \l\ shows the average spectrum of 
Mrk 335, in which the chosen continuum spectral windows 
have been highlighted. 

Two observations showing clearly spurious fluxes were re- 
moved; one in NGC 3516 (Julian day: 2449760.78), and one 
in NGC 5548 (Julian day: 2448993.89). No correction for red- 
dening was applied, as this should have no qualitative influence 
on our results (See Sect. 15. 3> . 
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Table 1. The 15 objects of the sample. The redshifts are taken from NED. The columns "SWP" and "LWP/R" give the number 
of observations in each wavelength range. The column "A!T" gives the monitoring duration in the rest frame (at 1300 A), while 
the average luminosity Logo in the band 1300-1350 A is given in the seventh column. 



Name 


Classification 


z 


SWP 


LWP/R 




Lnoo 


Eb-v 












(year) 


(xl0 40 ergs^) 




Mrk 335 


Q Pv f Pr t 1 


0.0257 


26 


28 


12.66 


622 58 +22 78 


0.059 


Mrk 509 


Spvfert 1 


0.0343 


39 


32 


15.06 


1 571 8 +67 16 


0.060 


Mrk 926 




0.0475 


22 


16 


14.51 


1 1 C2 Q +H7 1~\ 


0.053 


Mrk 1095 

1V±1 IS. 1U7J 


Spvfprt 1 


0.0331 


35 


23 


10.82 


Qf)S s 1 +79 44 


0.170 


NGC 3516 


Spvfert 1 


0088 


71 


22 


16.67 


40.51 ±2.44 


0.054 


NGC 3783 


Seyfert 1 


0.0097 


95 


84 


13.54 


70.01 ±2.13 


0.141 


NGC 4151 


Seyfert 1 


0.0033 


153 


137 


18.09 


26.45 ±2.17 


0.031 


NGC 4593 


Seyfert 1 


0.0083 


20 


15 


8.31 


12.65 ±1.22 


0.034 


NGC 5548 


Seyfert 1 


0.0171 


175 


148 


16.60 


189.70 ±5.21 


0.024 


NGC 7469 


Seyfert 1 


0.0163 


65 


15 


17.79 


202.78 ±4.59 


0.079 


3C 120.0 


Seyfert 1 


0.0330 


43 


21 


15.40 


159.96 ±9.47 


0.160 


3C 273 


Quasar 


0.1583 


124 


114 


15.36 


91782±1610.8 


0.027 


3C 390.3 


BLRG 


0.0561 


99 


11 


16.37 


292.98 ±16.11 


0.071 


Fairall 9 


Seyfert 1 


0.0461 


139 


63 


15.80 


1981.8 ±110.48 


0.042 


ESQ 141-55 


Seyfert 1 


0.0371 


26 


16 


11.33 


1785.1 ±101.40 


0.075 



The light curves at 1300-1350 A for all objects in our sam- 
ple are presented in Fig. [2 The monitoring durations in the rest 
frame of the objects are between 8.31 years (NGC 4593) and 
18.09 years (NGC 4151), see Table[T] 

4. Temporal analysis 

4.1. Structure function analysis of the light curves 

The first-order structure function of a light curve x(t) is a func- 
tion of the time lag r, and is defined by: 

SF,(T) = <(x(f + T)-x(f)) 2 >, (6) 

where (y) denotes the average of y over t. The structure func- 
tion (hereafter SF) an alysis was introduced in astronomy by 
ISimonetti et al] (|l985l). and is related to power density spec- 
trum analysis (Paltani 1999). It measures the amount of vari- 
ability present at a given time scale t. It has the advantage 
of working in the time domain, making the method less sen- 
sitive to windowing and alias problems than Fourier analysis. 
We shall use here the property that, if a maximum characteristic 
time scale T max is present in the light curve, the SF is constant 
for r > r max , reaching a value equal to twice the variance of 
x(t). Below T max the SF rises with a logarithmic slope of two at 
most. On very short time scales, the SF is dominated by the un- 
certainties on the light curve, and presents a plateau at a value 
equal to twice the average squared measurement uncertainty. 

We estimate the SFs of our light curves by averaging flux 
differences over predefined time bins, considering only the 
bins containing at least six pairs. We oversample the SFs (i.e. 
the bin-to-bin interval is smaller than the bin width) in order 
to emphasize their characteristics. The bins are geometrically 
spaced, i.e. the bin-to-bin interval is constant on a Log scale. 
Finally, we do not attempt to determine error bars on the SFs, 



as none of the prescriptions found in the literature see m sat- 
isfacto ry. For example, in the prescriptions of ISimonetti et alJ 
( 1985) and CP01, the uncertainty on the SF values is propor- 
tional to n7^ 2 , where n, is the number of pairs in bin i. These 
prescriptions produce underestimated error bars at large t (il- 
lustrated in Fig. 3 of CP01), because the number of pairs is 
increasing roughly exponentially with r because of the geo- 
metric spacing of the bins, while the total information in the 
light curve is finite. 

Structure functions of the 1300-1350 A light curves for the 
15 objects of our sample are presented in Fig. [3] The asymp- 
totic values at twice the light curve variances are also shown. 
Narrow structures in the SF (see e.g. Fairall 9, 3C 390.3) are 
very probably due to the high inter-correlation of the SF (be- 
cause the same measurements are used in several bins), and are 
probably not physical. 

In Fig. we observe that, for a majority of objects, the SFs 
show a plateau at about twice the variance of the light curves. 

However, in some cases, the SF continues to increase for t 

I — — • 1 J 71 

larger than ~ 5 years. Paltani et al. ( 1998) presented a similar 

analysis on 3C 273, and concluded that a second component, 
varying on long time scales, was present, particularly at long 
wavelengths. Such a component appears in the form of a SF 
rising sharply at large r. To cope with the possibility of the ex- 
istence of a secon d component, w e fit the SFs using the same 
function as in 

SF r (r) = 2e 2 + { A(t/7 "™ x) °'' t < T ™ x ) + Bl P , ( 7 ) 

[ A, T > T max J 

where e is the average uncertainty on the light curves, r max the 
maximum variability time scale measured at the start of the up- 
per plateau. A and a are the parameters of a first component of 
variability. B and B are the parameters of a second component. 
As we assume that this second component is slowly varying, 



4 



P. Favre et al.: AGN variability time scales and the discrete-event model 



GO 



6 - 



I 

■ i I 



3C 273 



I .... I .... I 



NGC 7469 







3C 390.3 






NGC 3783 



J 



3C 120.0 




i .... i 



60 




Mrk 335 




NGC 4593 




ESQ 141-55 




I .... I .... I .... I 



Mrk 926 








NGC 3516 




_1_ 



1980 1985 1990 1995 1980 1985 1990 1995 1980 1985 1990 1995 I960 1985 1990 1995 

Time (year) 



Fig. 2. Light curves in the range 1300-1350 A for the 15 objects of our sample. Only one error bar is drawn on each panel for 
clarity (in the upper- or lower-left corner). 



we adopt /3 — 2, the maximum value for the slope of an SF, 
which indicates very slow variability. 

We first fit the 1300-1350 A SFs; we shall discuss the 
longer wavelength SFs in Sect. 14.31 The values of T max are given 
in Table|2 corrected for time dilatation. The best fits are shown 
on Fig.[5]by a continuous line, e was fitted as a free parameter. 
The best fit values found for these parameters were compara- 
ble to the noise variance mesured in the light curves. In four 
objects (NGC 3516, NGC 4151, NGC 7469, Fairall 9), the fit 
does not converge unless we fix the noise parameter e and the 
slope a in the fits. We have tested that reasonable choices of 
e and a have little influence on the value of T max in those four 
SFs. We checked as well that the monitoring duration (AT in 
TableHJ had no influence on the measure of r max . 



The SFs of NGC 3516, and NGC 4151 show strong struc- 
tures between 0.1 and 1 year, but their behavior at r > 1 is 
compatible with an extrapolation of their behavior at r < 0.1. 
Our fit interprets the strong structures as evidence of a T max , 
but we need to check that the structures themselves do not re- 
sult from the sampling. To do that, we simulated light curves 
using a random walk, and projected them on the original light 
curve sampling. None of our simulations reproduced the ob- 
served structures in the SFs, and we consider therefore that a 
maximum time scale is really present in these two objects. The 
shape of the SFs forces us however to make the error bars on 
T m ax extend up to 1 year in these two objects. In Fairall 9, no 
plateau can be seen, but there is a clear change in the slope of 
the SF between r = 0.1 year and r = 3 years. We interpret this 
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Lag t (year) 



Fig. 3. Structure functions of the light curves at 1300-1350 A. The horizontal dashed lines show the asymptotic value of twice 
the variance of x(t). The continuous line shows the best fit, while the value of r max is shown with its uncertainty above the SF to 
ease readability. The vertical dot-long dash lines indi cate the location of the characteristic time scales found in CP01, while the 
short dash-long dash line indicates the value found in Paltan i et alJ (fl998>. 



as the presence of a r max in this range of time scale, followed 
by a very strong slowly varying component. In NGC 3783, a 
similar change of slope occurs, and it is impossible to locate 
Tmax unequivocally. This is reflected in the error bar on r max for 
this object. Several other objects show the existence of a second 
variability component, but it does not affect the measurement 
of r max . 

We repeated the analysis of the SFs without including a 
slow component of variability in the model, i.e. we fitted the 
data with B = 0. The values of r max found are compared in 
Table |3 The values found with B = are all inside the er- 



ror bars except for NGC 3516, NGC 3783 and NGC 4151. For 
NGC 3516 and NGC 4151, the fits do not represent the data. 
Thus for a majority of objects, the addition of a second com- 
ponent has no effect while it significantly improves the fits for 
NGC 35 1 6, NGC 3783 and NGC 4 1 5 1 . 

We conclude that with these data, we cannot decide if a 
second component is detected or not. This second variability 
component, while interesting per se, is outside the scope of this 
paper, and shall not be discussed further. 

For all the objects, we estimate the effect of the binning on 
T m ax by computing and fitting 100 SFs (at 1300-1350 A) for 
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Table 2. Maximum variability time scales r max of the 1300- 
1350 A light curves, corrected for time dilatation. The third 
column shows the values of T max obtained with B = 0, i.e. the 
SFs have been fitted without a slow variability component. The 
last column presents the corresponding range of event durations 
2yUi3oo, deduced from the simulations. 



Object 




1"max,B=0 


2/^1300 




(year) 


(year) 


(year) 



Mrk 335 
Mrk509 
Mrk 926 
Mrk 1095 
NGC 3516 
NGC 3783 
NGC 4151 
NGC 4593 
NGC 5548 
NGC 7469 
3C 120.0 
3C 273 
3C 390.3 
Fairall 9 
ESQ 141-55 



0.260±0.120 
0.550±0.100 
0.200±0.120 
0.400±0.140 
046+ 1 009 

vj.utu „ 06 
f) T, 1 Q+0.500 
U,J1 -031 

037 +1003 

"'" J -0.007 

0.284±0.230 

0.150±0.030 

0.022±0.005 

0.234±0.070 

0.452±0.050 

0.631±0.070 
341+2.717 

0.997±0.130 



0.222 
0.550 
0.227 
0.438 
4.896 
1.229 
2.200 
0.262 
0.154 
0.023 
0.287 
0.459 
0.635 
0.760 
1.014 



0.238±0.169 
0.489±0.381 

0.180+0.23?, 

0.242^ 
050 +L35 ° 
300 +L500 

_ 02g0 

040 +L50 ° 
0.686±0.452 
0.196±0.063 
0.035±0.014 
0.321+0.119 
0.559±0.117 
0.498±0.240 

350+ 2 - 800 

u.jju _ 0250 

1.630±0.816 



150 



B 

3 




100 



0.001 



0.01 0.1 1 

Input 2/u slm 1300 value (year) 



Fig. 4. Histogram of the event durations 2/4^1300 input to the 
simulation that produce a measured T max in the range 0.4-0.5 
years, for 3C 273. The inset shows the same diagram in linear 
scale, fitted with a Gaussian. 



each object with different binnings, with corresponding bin-to- 
bin intervals between 0.002 year and 0.1 year. For all objects, 
the distribution of the measured T max is mono-peaked, meaning 
that a single value of T max was always found by the algorithm. 
The width of this peak determines an empirical uncertainty on 
r max (Table |3. 

We measure r max from 0.022 to 0.997 year for the 15 
objects of our sample. Our time scales are of the same 
order of magnitude as the one found by previous stud- 
ies in the optical-UV JHooketalJ fl994 j iTrevese etalj|l994[ 
ICristiani et ail 1 199(4 IPaltani et alJ Il998t Idiveon et all Il999l 
CSV0 0). The time scales found by CP01 and IPaltani et alJ 
dl998b are indicated in Fig.[3]by a vertical line and are in rea- 
sonable agreement with what we have found, 3C 390.3 and 
NGC 3783 excepted. The discrepancy can be explained by the 
fact that CP01 "detrend" their SFs with a linear component (i.e. 
they remove a linear fit from their light curves), arguing that 
the measured SF will deviate from their theoretical shape if the 
light curve shows a linear trend. Our method of including a 
second, slowly variable component in the structure functions is 
more general than the "detrending", because it makes less strict 
assumptions on the temporal properties of the slowly varying 
component. It is nevertheless equivalent in the case where a lin- 
ear trend is effectively present in the data. Our method is also 
more consistent in the sense that all components are handled 
in a similar way. Furthermore, a linear trend would make little 
sense in objects like NGC 4151. 



4.2. Relation with the event duration 

An easy way to explain the existence of a maximum variabil- 
ity time scale is provided by the discrete-event model. For a 
Poissonian seque nce of events, the SF is proportiona l to the SF 
of a single event (Palta nJl996tlAretxaga et all 19971 CSV00), 
and only has structures on time scales shorter than the event 
duration. We interpret the observed SFs using discrete events 
for which we assume a triangular, sym metric shape. Th e event 
shape was chosen for its simplicity ( Palta ni et al Jl 19981) . but it 
has been shown that choosing other shapes does not affect sig- 
nificantly the results given by the temporal analysis (CSV00). 
Our events are described with only two parameters at wave- 
length A: the event amplitude Hx and the event duration 2fiy, 
being defined as the time needed to reach the maximum flux. It 
follows that £l = 1/2 and k\ = 1/3 in this case (see Sect.|2j. 

While the SFs are in theory able to measure the event du- 
ration, this measure can be affected by the noise and sampling 
of the light curve in a complex and unpredictable way. To test 
if T"max measures a property of the light curves, and not of the 
sampling, we produce synthetic light curves by simulation, and 
measure their r max . In the simulations, we add randomly trian- 
gular events with a given duration 2 / u simI3 oo, keeping the same 
sampling as the original light curve. A noise with an amplitude 
equal to the average of the instrumental noise is added to each 
light curve. We take 40 test values for 2yU s i m ,i3oo, from 10~ 3 to 
10 years, geometrically spaced. 

For each object and each 2/z s i m> i3oo, we build 1000 light 
curves, compute their SFs and, measure T max using the method 
described above. The event rate N is randomly chosen between 
5 and 500 events per year. The result of the simulation is, for 
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(e) 2100-2150 A 




(d) 1850-2000 A 




(c) 1700-1750 A 




0.1 1 

T(years) 



o.oi 
l 




(1) 2975-3025 A 





(j) 2700-2750 A 




(i) 2550-2600 A 




(h) 2425-2475 A 



^g£j 




0.1 1 

T(years) 



Fig. 5. SFs for the light curves 1300 to 3000 A, for NGC 5548. The horizontal dashed lines show the asymptotic value of twice 
the variance of the light curve. The best fits are shown by a continuous line. 



each object and each 2/z s i m i3oo, the distribution of the result- 
ing r max . These distributions allow us to determine which input 
2/i S im,i300 can provide the observed r max . For 3C 273 for exam- 
ple, the distribution of 2 / u s i m i3oo produced a peak around 0.56 
year, as represented on Fig.|4] The fit of the peak of Fig.|4]with 
a Gaussian gives 2jU s j m; i3oo = 0.559 + 0.1 17 years. 

We note that we never observe a B parameter (see Sect. 14. It 
significantly larger than in our simulations. This is expected 
as we do not include a second component. 

The simulations for all objects showed a result qualitatively 
identical to that for 3C 273, i.e. the distributions of 2// s ; m i3oo 
present a single peak. This means therefore that, for each ob- 
ject, r max determines a unique event duration, that can be de- 



rived from the simulations. The values of 2// 1300 are given in 
Table [2] We stress that our simulations are driven by the real 
sampling of the light curves, and ar e there fore more specific 
than, for example, those discussed in IWelshl 179991 or CP01. 



4.3. Wavelength dependence of the event duration 

For each object of the sample, we apply the method described 
in Sect. l4.ll to compute the event durations from the light curves 
1450 to 2975 A. 

Structure functions from 1300 A to 3000 A are presented 
in appendix A for each object, along with a description of the 
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Fig. 6. r max as a function of wavelength for the 15 objects. The T max are corrected for time dilatation. The two objects which 
present a very strong increase in T max at long wavelengths are shown with a dashed line. 



particularity of each set of SFs. As an example, we show the 
case ofNGC 5548 in Fig. 

It is not possible to deduce a value of T max for all 180 light 
curves. This is mainly due to the fact that, for some of the ob- 
jects, the number of observations in the LW range is too small. 
In addition some particular light curves are very noisy. In such 
cases, the SF usually does not have the canonical shape, and the 
fit does not succeed. We thus reject the time scale correspond- 
ing to those particular SFs. Each individual case is described 
in appendix A. In some cases, the noise is such that, although 
Tmax can be derived from the SF, its uncertainty is very large. 
In such cases, one should interpret any variation in r max with 
caution. 



Fig- ED presents the variability time scale r max as a func- 
tion of the wavelength for all the objects. We find that r max is 
reasonably constant over the wavelength range we use, as the 
small fluctuations can be explained by the difficulty to measure 
a precise value of r max on some noisy SF. For the particular 
case of NGC 5548 for example, the variations are inside the 
uncertainties d erived in the previo us section. This result was a 
also found by Palt ani et al.l dl998h for 3C 273. 

In two objects however, NGC 7469 and Fairall 9, the SFs 
present a very strong increase in r max at long wavelengths. For 
NGC 7469, this is due to a lack of short term sampling of the 
LW light curves, which prevents the recovery of any time scale 
below 0.5 year. In Fairall 9, a similar lack of short term sam- 
pling affects the determination of r max . However, the values of 
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10 45 



Fig. 7. Variability time scales 2yUi3oo as a function of the lu- 
minosity of the objects in the range 1300-1350 A. The BCES 
regression is shown with the dashed line. 



Tmax m me LW range are within the uncertainties on r n 
mined at 1300 A. 



deter- 



5. Interpretation in terms of discrete-event model 

5.1. 2//1300 as a function of the luminosity 

The relationship between the event duration and the luminosity 
has some important consequences for the discrete-event model 
that we will discuss below (Sect. \5A\ . We note, however, that 
only a weak dependence between the event duration and the 
average luminosity of the object is possible as the values of the 
former cover less than two orders of magnitude while the latter 
covers four orders of magnitude. We thus measure a physical 
time scale which has at most a small dependence on the lumi- 
nosity of the objects. 

The event duration 2yUi3oo as a function of the average lu- 
minosity Li3oo of the object at 1300-1350 A is shown in Fig.0 
We use a simple cosmology with Ho = 60 km s -1 Mpc~\ and 
qo = 0.5 throughout. 

Using Spearman's correlation coefficient, we find a correla- 
tion between the event duration and the luminosity (s = 0.38), 
marginally significant at the 16% level (Null hypothesis). The 
dependence of 2//1300 on Looo can be expressed as 2yifi3oo oc 
^Boo , where the index 6 - 0.21 ±0.11 has been de termined 
using the BCES linear regressions JAkritas & Bershad\ll 996). 



5.2. The steady component C A 

The steady component Ca (Sect. [3 can have various physical 
origins. For example, it can be associated to the non-flaring 



part of the accreti on disk, or to the host-galaxy stellar con- 
tribution (|cidFernandeIetaO |l993), CSV00). We shall how- 
ever continue our discussion in a model-independent way. We 
can constrain Ca for a particular light curve by imposing that 
it does not exceed the minimum observed luminosity i™ n . 
On the other hand, Ca — i s an obvious lower limit (al- 
though Paltani & Walterj dl996l) argued that C A > 0, at least 
for A > 2000 A). In the following, we shall use these two con- 
straints as limiting cases. 

5.3. Spectral shape of the event amplitude, event 
energy and event rate 



For each object, we compute the event amplitude Ha from 
Eq. J5} using both Ca =0 and Ca = L™ m . Fig.|S]shows the spec- 
tral shape of Ha for NGC 5548. We integrate E A = 2/j A k L HA 
interpolated over the wavelength range 1300-3000 A to obtain 
the energy E released in one event, assuming isotropic emis- 
sion: 



,-3000 A 
Jl300 A 



-3000 A 

E = I /.^d. I . 



(8) 



We use the event duration at 1300-1350 A, as it can be consid- 
ered constant over the wavelength range considered. The event 
rates N are derived from Eq. ©. 

Table |3 gives the event energies and rates found with the 
Ca = 0, and Ca = L™ m assumptions. Event energies are found 
in the range 10 4S - 10 52 erg, and maximum event rates in the 
range 9-1 133 event year -1 (Ca = 0) while minimum event rates 
are found in the range 2-270 event year -1 (upper limit of Ca). 
Fig.|5]shows the event energy £ as a function of the 1300-1350 
A luminosity, using the upper and lower limits on Ca- 

With Ca — 0, E and L1300 are clearly correlated (Spearman 
correlation coefficient: s = 0.85), with a Null hypothe- 
sis probability of less than 0.1%. A linear regression gives 
E oc Li3oo 7 , with y = 9 6 ± 0.09 using the BCES method 
dAkritas & Bersha dvl 19961) (Fig.|5] left panel). 

When using Ca = the correlation is preserved 

(Spearman: s = 0.89), with a Null hypothesis probability of 
less than 0.1% and an index y — 1 .02+0. 1 1 (Fig. [9] right panel). 
We note that E and Looo are lower limits since no corrections 
for reddening were applied. The correlation should be pre- 
served by applying the corrections, as both E and Z4300 would 
be scaled by the same factor. We tested this hypothesis by cor- 
recting the light curves for reddening and recomputing the re- 
lation E = /(L1300). The correlation is preserved (Spearman 
correlation coefficient: s = 0.83) with an index y = 0.93 + 0.10 
for the case Ca = as well as for the case Ca = £™ n (Spearman 
correlation coefficient: s = 0.84) with an index y = 0.97 ±0.11. 
These values are consistent with the non-dereddened values. 

The event rates as a function of the object luminosity 
are given in Fig. [K)] for both the upper and lower limits on 
Ca- A non-significant anticorrelation is found in the Ca — 
case (Spearman's s = -0.12, Null hypothesis probability of 
66.64%), while a marginally significant anticorrelation is found 
for Ca = £™ n (Spearman's s = -0.39, Null hypothesis 
probability of 15%). Using the data corrected for reddening, 
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Table 3. Event rate N and energy E for each object, with the Ca = 0, and Ca = L™ m assumptions. 



Name 


N 




E 


N 


E 






(year 




(xlO 50 erg) 


(year -1 ) 


(xlO 50 


erg) 






C A = 





C A = 


L min 

,( 




Mrk 335 


160.76 ± 


123.78 


0.25 ± 0.26 


17.82 ± 41.22 


0.72 ± 


0.76 


Mrk 509 


39.28 ± 


31.89 


2.03 ± 0.67 


9.09 ± 15.34 


5.08 ± 


1.64 


Mrk 926 


24.84 ± 


17.20 


3.94 ± 0.87 


14.66 ± 13.22 


6.12 ± 


1.30 


Mrk 1095 


149.32 ± 


90.68 


0.26 ± 0.08 


18.45 ± 31.88 


0.99 ± 


0.35 


NGC 3516 


104.90 ± 


23.42 


0.024 ± 0.002 


75.44 ± 19.86 


0.037 ± 


0.004 


NGC 3783 


50.11 ± 


48.48 


0.08 ± 0.02 


28.21 ± 36.37 


0.11 ± 


0.03 


NGC 4151 


32.12 ± 


8.62 


0.05 ± 0.002 


31.16 ± 8.49 


0.056 ± 


0.002 


NGC 4593 


10.32 ± 


7.09 


0.08 ± 0.02 


2.32 ±3.36 


0.19 ± 


0.05 


NGC 5548 


51.59 ± 


18.49 


0.20 ± 0.03 


29.52 ± 13.99 


0.30 ± 


0.05 


NGC 7469 


1133.8 ± 


543.73 


0.02 ± 0.01 


269.88 ± 265.39 


0.04 ± 


0.02 


3C 120.0 


27.53 ± 


12.42 


0.45 ± 0.48 


13.22 ± 8.61 


0.65 ± 


0.67 


3C 273 


62.44 ± 


14.51 


71.20± 27.34 


10.24 ± 5.88 


189.62 i 


53.83 


3C 390.3 


8.94 ± 


4.82 


1.64 ±0.49 


5.86 ± 3.90 


2.49 ± 


0.78 


Fairall 9 


8.82 ± 


6.44 


11.11 ± 1.49 


6.50 ±5.53 


13.03 ± 


1.74 


ESQ 141-55 


9.75 ± 


5.11 


9.38 ± 1.85 


2.69 ± 2.69 


20.03 ± 


3.86 
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Fig. 9. Event energy as a function of the object luminosity in 
the hypothesis Ca — (left panel), and using the upper limit 
of Ca (right panel). The BCES regression is shown with the 
dashed line. 



Fig. 8. Event amplitude spectrum H\ as a function of wave- 
length for NGC 5548, with C A = (squares), and C A = Lf n 
(open triangles). The lower curve has been slightly shifted to 
the right to ease readability. 



one finds as well a non-significant anticorrelation (Spearman's 
s — -0.02, Null hypothesis probability of non-correlation of 
93.96%) in the case Ca = 0, and a non-significant anticorre- 
lation for Ca = £™ n (Spearman's s — -0.17, Null hypothesis 
probability of 54.12%). 



5.4. Constraints on the variability-luminosity relation 



IPaltani & Courvoisierl i 19941) showed that the variability of 
a similar sample was anticorrelated to the object luminosity. 
PC97 confirmed this result in the rest-frame of the objects and 
found crf^ ~ with 77 = -0.08 ± 0.16. 

We showed in Sect. 15. II that the event duration 2yUi3oo is 

a shallow function of the luminosity and in Sect. 15.31 that the 

y 

measured event energy E oc Z4300 , with y 1. These two 
results, expectedly, lead to values of rj in agreement with the 
measure of PC97. We thus established that the event parameter 
which drives the cr(L) dependence is the event energy, and not 
its duration, nor its rate. 
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Fig. 10. Event rate as a function of the object luminosity in the 
hypothesis Ca = (left panel), and using the upper limit of Ca 
(right panel). 

Table 4. Slope r\ of the variability-luminosity relation for the 
two C,\ assumptions (see text). 







r 




6 


n 




Ca = 
C A = LT 


0.96 
1.02 


±0.09 
±0.11 


0.21 
0.21 


±0.11 
±0.11 


-0.13 ± 
-0.10± 


0.06 
0.05 



PC97 showed that in this case one should expect a 
variability-luminosity relation in the form: 



o-(L) oc L 130 o 



(9) 



We present in Table 0]the different values of the slope 77 of the 
variability-luminosity relationship which are consistent with 
the value ofPC97. 



6. Discussion 

6.1. Event duration and black-hole physical time 
scales 

We have shown that a characteristic variability time scale ex- 
ists, which can be measured in the light curves. It can be as- 
sociated with the event duration in a model-independent way. 
We have obtained event durations in the range 0.03 to 1.6 
years, which may possibly be related to the four physical 
time scales associated to black holes. They all depend on the 
black hole mass Mbh and Schwarzschild radius We re- 
view them below, from the f astest to the slowest, following 
lEdelson & Nandral d 1 999t) and lManmoto et all dl99fl 

1. The light crossing time is given by t\ c = 3.01 x 
10- 5 M 7 (fl/10i? s )year; 

2. the ADAF accretion time scale (comparable to free-fall ve- 
locity) is given by f acc > 4.38 x 10- 3 M 7 (#/100fl s ) 3/2 year; 

3. the gas orbital time scale is given by f or b = 9.03 x 
10- 4 M 7 (fl/10fl s ) 3/2 year; 




0.001 



0.1 

Time (year) 



Fig. 11. Time scales vs black hole mass. The gray area shows 
the range of timescales found in this study. 



4. the accretion disk thermal time scale f t h = 1.45 x 
10- 2 (0.01/a visc )M 7 CR/10fl s ) 3/2 year, 

where R is the emission distance from the center of mass, 
M 7 = M/1O 7 M , and the Schwarzschild radius is defined 
Rs = 2GM/c 2 . f t h and f or b can produce the range of time scales 
observed here for reasonable black hole masses (see Fig. llU . 
However, all these time scales depend linearly on the mass of 
the object, hence on the luminosity. The lack of strong depen- 
dence of the event duration on the object luminosity allows us 
to exclude all these mechanisms as likely candidates for the 
origin of the variability. 

6.2. Physical nature of the events 
6.2.1 . Supernovae 

In the starburst model dAretxaga et alJll997l) . the variability of 
AGN is produced by supernovae (SNe) explosions and compact 
supernovae remnants (CSNRs). The SNe generate the CSNR 
in the inter action of their ejecta w ith the stellar wind from the 
progenitor. lTerlevich et al J 11992) showed that the properties of 
CSNR match the properties of the broad -line region of AGN. 

lAretxaga & Terlevichl dl993l 1 1994 modeled the B band 
variability of the Seyfert galaxies NGC 4151 and NGC 5548 
with this model. For NGC 4151, an event rate of 0.2-0.3 events 
year 1 was found. However, typical prediction s of the m odel 
are m ore of the order of 3-200 events year -1 dAretxaga et alJ 
119971) . consistent with what we found. The event energy has to 
be constant (3-5xl0 51 erg; e.g. lAretxaga et al.1dl997h ). in clear 
contradiction with our result. The lifetime of CSNR (0.2-3.8 
years) is compatible with the event durations found here. But 
no correlation with the object luminosity is expected, the more 
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luminous objects simply having higher SN e rates . Again, this 
is contrary to our results. Finally, Aretxag aetalJJl997t) show 
that a -1/2 slope of the <x(L) relation should always be found 
with this model, which again is not observed. 

6.2.2. Magnetic blobs above an accretion disk 

In this model, each event is associated with the discharge of an 
active magnetic blob a bove an accretion disk. In the model pro- 
posed bv lHaardt et al.l i 19941) . a fraction of the local accretion 
power goes into magnetic field structures allowing the forma- 
tion of active blobs above the disk. Reconnection of the mag- 
netic field lines in the corona permits the transfer of the energy 
into kinetic energy of fast particles. The energy is stored and 
released in the so-called charge and disc harge times t c and U 
with ?d <sc t c . Using th e dynamo model of [ Galeev et alJ (fl979) 
for the blob formation, Haar dt et alJ (11994) show that f<j scales 
with the blob size which itself scales with the total luminos- 
ity L of the source. This trend is clearly not seen in our data 
(Fig-ID- Furthermore, the total number of active loops iV tot , at 
any time, does not depend on the luminosity nor on the mass of 
the object. The blob rate iVtotAd becomes therefore proportional 
to Zr 1 , also in clear contradiction with our results that show that 
the event rate is not correlated with the luminosity. Finally, the 
energy released by a single blob can be written E - fd^biob, 
where Lbiob is given by Eq. (7) of H aardt et alJ 11994). The en- 
ergy E released then goes with L 2 , also in contradiction with 
the results deduced here in which E oc L. 



6.2.3. Stellar collisions 

ICourvoisier et alJ dl996l) proposed that the energy radiated in 
AGN originates in a number of collisions between stars that 
orbit the supermassive black hole at very high velocities in a 
volume of some 100 Rs,. They computed the rate dn/dt of head- 
on stellar collisions in a spherical shell of width dr, located at 
distance r from the central black hole. The stars are assumed to 
have mass Mq and radius R of the Sun. This reads: 



dn o 

— = (p 2 4nr 2 dr)v K 7TR 2 

dt 



(10) 



where p repr esents the ste ll ar den s ity and vk the K eplerian 
veloci ty dCourvoisier et al.l dl996l) : iTorricelli-Ciamponi et alJ 
(2000)). Assuming that the stars are located in a shell of in- 
ner radius a and outer radius b and are distributed following a 
density law p = No(r/ro)~ a ^ 2 , with slope a where No and ro are 
constants, we find for the kinetic energy released by one event: 



(id 



f d f £drv K (r)r 2 -« 

where a, a and b are parameters of the stellar cluster. Using 
v K = y/GM tot /r leads to: 



E = 



\M Q G$ b a drMl t A 



J a 



(12) 



Neglecting the cluster's mass with respect to the black hole 
mass, i.e. making the assumption that M tot (r) ~ Mbh, we fi- 
nally have: 



E = ^M Q GM BH f(a,a,b) 



(13) 



where / is a function of the cluster parameters only, indepen- 
dent of Mbh- 

We need now to relate the average luminosity to the black 
hole mass. This relation comes from L — J \ jlMQv\dnldt = 
(M G R 2 Q 7T 2 G V2 N r%)M% 2 f(a,a,b). As £ oc M BH and I oc 
Mg{?, we finally have: 



E oc V 



(14) 



drM? ot r^~ 



This relation is a relatively good approximation of the trends 
seen in Fig.|5] although not completely satisfactory. 

In this model, the variability time scale 2/^ is expected to 
be related to the time needed to the expan ding sphere to become 
optica lly thin. This point is discussed in ICourvoisier & Ttirlerl 
(2004) who found that for clumps of about one Solar mass, 
the expansion time is about 2 x 10 6 seconds. This enters in the 
range of time scales found here. The collision rate should be 
going with M^ 2 , which implies N oc l' /3 , which is not seen in 
our data. 



6.2.4. Other models 

For the sake of completeness, we note that gravitational micro- 
lensing models, in which populations of planetary mass com- 
pact bodies randomly cross the line of sight of an observed 
AGN, can be invoked to explain the long term variations over 
several years. But in the case of low-redshift Seyfert galaxies, 
which forms the majority of ou r sample, the p robabilities of 
micro-lensing are not significant dHawkinsl200ll) . 

Final ly, a model of accretio n disk instabilities has been sug- 
gested ( Kawaguch i et al . 1998) to explain the optical variabil- 
ity of AGN. The SOC state model IMineshige et al.lll994l) is 
producing power density spectra in good agreement with the 
observations but, since it is not an event-based model, it is dif- 
ficult to use the measurements discussed here to constrain it. 



7. Conclusion 

We showed the existence of a maximum variability time scale 
in the ultraviolet light curves of 15 Type 1 AGN, in the range 
1300-3000 A. We found variability time scales in the range 
0.02-1.00 year. 

In the framework of the discrete-event model, we showed 
that these time scales can be related to the event duration in a 
simple manner. A weak dependence of the event duration with 
the object luminosity at 1300 A is found. The event duration is 
not a function of the wavelength in the range 1300 to 3000 A. 

The event energy per object varies from 10 48 to 10 52 erg 
with a corresponding event rate comprised between 2 and 270 
events per year, assuming the presence of a constant component 
in the light curves. 
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Our results do not depend on the constant component Cx- 
While we can only provide lower and upper bounds on Cx, its 
choice does not change the conclusions. 

The event energy is strongly correlated with the object lu- 
minosity. We show that the combined relations of the event en- 

1.02 0.21 

ergy E oc Z4300 , and event duration 2/ii3oo °c £4300 with 
the object luminosity, lead to the trend seen in the variability- 
luminosity relationship in the rest frame, i.e. that both variables 
are correlated with a slope of about 0.08. We thus established 
that the event parameter which drives the cr(L) dependence is 
the event energy, and not its duration, nor its rate. 

These results allow us to constrain the physical nature of the 
events. We show that neither the starburst model nor the mag- 
netic blob model can satisfy these requirements. On the other 
hand, stellar collision models in which the average properties 
of the collisions depend on the mass of the central black hole 
may be favored, although the model will need to be improved 
as the result we found (for instance, the lack of correlation be- 
tween the event rate and the luminosity) does not match the 
predictions. 

Acknowledgements. SP acknowledges a grant from the Swiss 
National Science Foundation. 
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